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Abstract 

A fluid flow and heat transfer model has been developed for the reactive, porous bed of the biomass gasifier to simulate pressure drop, 
temperature profile in the bed and flow rates. The conservation equations, momentum equation and energy equation are used to describe 
fluid and heat transport in porous gasifier bed. The model accounted for drag at wall, and the effect of radial as well as axial variation in 
bed porosity to predict pressure drop in bed. Heat transfer has been modeled using effective thermal conductivity approach. Model pre¬ 
dictions are validated against the experiments, while effective thermal conductivity values are tested qualitatively using models available 
in literature. Parametric analysis has been carried out to investigate the effect of various parameters on bed temperature profile and pres¬ 
sure drop through the gasifier. The temperature profile is found to be very sensitive to gas flow rate, and heat generation in oxidation 
zone, while high bed temperature, gas flow rate and the reduction in feedstock particle size are found to cause a marked increase in pres¬ 
sure drop through the gasifier. The temperatures of the down stream zones are more sensitive to any change in heat generation in the bed 
as compared to upstream zone. Author recommends that the size of preheating zone may be extended up to pyrolysis zone in order to 
enhance preheating of input air, while thermal insulation should not be less than 15 cm. 

© 2007 Elsevier Inc. All rights reserved. 
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1. Introduction 

The flow of fluids through porous media is of great prac¬ 
tical importance in many diverse applications, including 
the drying the wet biomass particles or grain seeds, the gas¬ 
ification of biomass or coal, ground water movement, 
regenerative heat exchange, surface catalysis of chemical 
reactions etc. In all these instances it is necessary to predict 
design parameters such as pressure drop, friction factor, 
heat transfer coefficients in order to predict optimum oper¬ 
ating conditions, size and design of equipment. In the past 
there has been a steady effort to improve the knowledge in 
this area. In most cases, Darcian flow assumed; however, 
some researchers have concentrated into non-Darcian flow 
as well. The pressure gradient across the porous bed is a 
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function of system geometry, medium porosity, permeabil¬ 
ity and physical properties of working medium. Consider¬ 
able progress has been made in past in establishing the 
relationship among velocity versus pressure drop or Rey¬ 
nolds number versus friction factor. 

Ergun (1952), Ranz (1952), Brownell et al. (1956), Tall- 
madge (1970), Comiti et al. (2000) developed dimensionless 
form to establish the relationship in frictional factor and 
Reynolds number to describe the fluid transport through 
porous media as 

./pore = a/ftepore + fi/R- e p 0 K + ^ (1) 

here, Re pore [ = pv 0 t/ p //(( 1 - £bed)] is the pore Reynolds 
number, / pore is the friction factor, a, [l and 1 are the con¬ 
stants. Except Tallmadge (1970), the value of /(is fixed at 0. 

Darcy law deviates as fluid velocity increases, since iner¬ 
tial contribution starts dominating. Thus in order to 
account the solutions for the complete regime of flow, 
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Nomenclature 



A 

area (m 2 ) 

£bed(T) 

radial bed porosity 

CV 

control volume 

T s 

sphericity 

dp 

current particle size (m) 

/l 

constant 

k 

thermal conductivity (W/m K) 

<7 

Stefan-Boltzmann constant 

m 

mass flow rate (kg/s) 

^bed 

Bed porosity (bulk) 

^Tuy 

number of tuyers 

a, /? 

constants 

2gen 

heat generation (kW) 

£ 

correction factor at reactor wall 

T 

temperature (K) 



Y 

mass fraction 

Subscripts!superscripts 

A m 

mean area (m 2 ) 

A 

ambient 

d 

particle diameter (m) 

bd 

dry biomass 

D 

geometrical diameters (m) 

char 

char 

K 

pressure drop parameter 

con 

contact of particles 

yyigm 

mean gas flow in each CV 

gen 

generation/released (heat) 

A P 

pressure drop (N/m 2 ) 

preheat preheating zone 

R 

gasifier tube radius (m) 

sg 

solid-gas-fluid contact 

T a 

ambient temperature (K) 

Tuy 

tuyer 

T 

thickness (m) 

oxd 

oxidation 

cu 

constants 

air 

air 

db 

initial particle diameter (m) 

dry 

drying zone 

A 

tube diameter (m) 

cl 

unit cell 

/ 

length (m) 

Dev 

developing flow 

A m s 

solid mass conversion 

pyr 

pyrolysis 

Q 

heat (kW) 

Pg 

producer gas 

Re 

Reynolds number 

s, g 

solid, gas 

AT 

temperature drop in unit cell 

vol 

volatile 

v 0 

superficial velocity (m/s) 

red 

reduction 



ash 

ash 

Greek letters 

eff 

effective 

P 

viscosity (kg m -1 s -1 ) 

ct 

contact of two particles 

^ex 

extinction coefficient 

Fd 

fully developed flow 

£bed 

average bed porosity 

pore 

pore 

p 

density (kg/m 3 ) 

rad 

radiation 

f 

friction factor 

Top 

open top 

0 

angle 

w 

water vapour 

(^ex) 

scaled extinction coefficient 




Forchheimer proposed the first expression in 1901, which is 
based on dimensional correlation (Churchill, 1988). 

— AP/l = a*pv 0 + AP/1) (la) 

where a* and /f are the constants. Ergun (1952) combined 
the Carman-Kozeny equation, and Blake-Plummer equa¬ 
tion in form of Forchheimer equation. Macdonald (Chur¬ 
chill, 1988; Kaviany, 1995) successfully used the 
Forchheimer’s equation and correlate data separately for 
different types of porous media in terms of various param¬ 
eters including particle surface roughness. In order to 
implement these models to the chemical reactors/gasifiers, 
uncertainties still exist- (i) varying tube-to-particle diameter 
ratios and (ii) drag due to reactor wall. Various authors ad¬ 
dress these issues for packed bed (Mehta and Hawely 
(Churchill, 1988); Cohen and Metzner (1981); Winterberg 


and Tsotsas (2000a); Winterberg and Tsotsas (2000b).) 
Mehta and Hawely describe the wall drag by taking the 
surface area of the tube into account in the definition of 
the hydraulic diameter, while Winterberg and Tsotsas 
(2000a,b) handle the radial porosity variation and wall 
drag in packed tubes. Sharma et al. (2004) developed an 
interesting fluid analysis to relate air and gas flow 
through the porous bed of biomass gasifier, with increasing 
fluid flow and decreasing particle size in the direction of 
flow. 

Extensive work has been carried out on heat transfer in 
packed bed (Singh and Kaviany, 1992; Cheng et ah, 1999; 
Kikuchi, 2001; Sharma, 2006). These include many 
approaches. The approach using effective thermal conduc¬ 
tivity, k eB , of the packed bed has been most popular in the 
recent literature (Chen and Churchill, 1963; Kaviany, 
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1995; Slavin et al., 2002; Kamiuto and Yee, 2005; 
Sharma, 2003, 2005). Slavin et al. (2002) have used the 
measurable parameter, surface roughness of the particle 
to avoid the use of any adjustable parameters in the deter¬ 
mination of k ef{ for uniform metallic spheroids hlled with 
ideal gas with no compression of particles. Sharma (2003) 
considered a packed bed of biomass and char filled with 
gas mixture, with surface contact between particles due 
to compression, Sharma (2005) further extended the Ar eff 
model for radiation transport by modifying shape factor 
using the extinction coefficient and also developed the 
heat transfer model using single control volume applied 
to the bed of gasifier. 

Logtenberg et al. (1999) have reported a fluid flow and 
heat transfer CFD model using 3-D Navier-Stocks equa¬ 
tions. Marafie and Vafai (2001) reported a non-thermal 
equilibrium, two-equation model to represent the fluid 
transport and energy transport. Jamialahmadi et al. 
(2005) presented a mathematical model for the prediction 
of pressure drop, bubble size, gas hold up and heat transfer 
accurately. 

In the present contribution, the focus is on analysis of 
heat transfer and fluid flow considering axial, radial varia¬ 
tion in bed porosity along with the effect of wall drag in 
order to predict accurate various flow rates, pressure drop, 
temperature profile through the gasifier bed and apportion¬ 
ment of air inflow from open top to tuyers. 


2. System description and experiments 

The open top downdraft biomass gasifier of capacity 
20 kWe developed in Indian Institute of Science, Bangalore 
(Sharan et ah, 1997) marketed by NETPRO has been cho¬ 
sen for the present work (Fig. 1). In this design, air enters 
from top as well as from radial tuyers and the producer gas 
exits from the bottom of the reactor. The gas is than made 


Biomass An- 


Thermal Insulation 


O Producer gas 

Preheating & 
Annular jacket 
zone 



Fig. 1. Sectional view of IISc-DASAG gasifier. 


to pass through an annulus around preheating zone so as 
to utilize a part of its heat for preheating of biomass. 
The outer wall of the gasifier is insulated to improve the 
thermal efficiency of the reactor. The system using pro¬ 
ducer gas for motive power generation, the gasifier should 
be coupled with a suitable cooling and cleaning unit. The 
blower has been used in the gas pipeline to create required 
suction to draw the gas. 

In order to validate the model, the experimental data 
has been generated for IISc-DASAG gasifier. The gasifier 
was ignited and time was allowed (about 20 min) for the 
initial transients to subside. Then the temperatures at var¬ 
ious locations inside the reactor were recorded using cali¬ 
brated K-type (Chromel-alumel) thermocouples. Two 
single-channel digital temperature indicators, a 12-channel 
selector switch were used in the present work to read out 
the temperature values at various locations. The pressure 
at tappings PI is recorded using U-tube manometer 
arrangements. Flow rate of producer gas was measured 
using a calibrated venturi meter. The procedure was 
repeated for different gas flow rates. The experimental 
set-up of a gasifier system has been shown schematically 
in Fig. 2. 

3. Formulation 

The packed bed in the gasifier has been modeled as a sat¬ 
urated porous medium in which fluid flow increases in the 
direction of flow due to gasification of the particles consti¬ 
tuting the bed. The flow of air and biomass consumption in 
a biomass gasifier is closely coupled with thermo-chemical 
processes. For simplicity, the gasifier is divided into six 
zones, each one corresponding to a dominant phenomenon 
as: preheating, drying, pyrolysis, oxidation and reduction 
zones. The annular jacket is considered as a separate zone. 
The preheating, annular jacket zone and drying zones are 
further divided into 3 control volumes (CVs) each; and 
the pyrolysis zone is divided into 2 CVs. Oxidation and 
reduction zones are constituted by one control volume 
each. 

3.1. Assumptions 

The porous bed is assumed to be isotropic, where solid 
and gas is considered to be in thermal equilibrium. Within 
each zone, the particle size, bed porosity and temperature 
are assumed to be uniform. The following assumptions 
have been invoked for each zone: 

Drying!preheating zone'. The particles are taken to be 
spherical with equal diameters. Complete removal of 
moisture is assumed in this zone. No shrinkage in parti¬ 
cle size. 

Pyrolysis zone : All volatiles are assumed to be released 
in this zone and moving into gaseous phase leaving 
behind char and ash in the solid phase. Heat of the pyro¬ 
lysis reactions is taken to be negligible. 

































A.Kr. Sharma / Int. J. Heat and Fluid Flow 28 (2007) 1518-1530 


1521 


Biomass Air (Top) 



Downdraft (Biomass) Spray Coolers Sand Bed filters Blower Valve Venturi Burner 

IISc-NETPRO Gasifier 


Fig. 2. Schematic diagram of a gasifier system. 


Oxidation zone'. Oxygen in available air is assumed to be 
completely consumed in this zone. Char oxidation is 
negligible in comparison with volatile oxidation. 
Reduction zone : All char is taken to be consumed in this 
zone. The ratio of heat generation in oxidation zone to 
heat absorption in reduction zone is fixed at 2.8. 

3.2. Mathematical equations 

3.2.1. Mass balance 

In the biomass gasifier, the consumption of biomass, 
water vapour, flow of air from open top and through the 
tuyers can be related with gas flow rate at the gasifier exit 
using global mass balance as 

fffur — 77?Top d UTuyf^Tuy — ^bd 7W aS h ( 2 ) 

During gasification process, the solid phase converts pro¬ 
gressively into gaseous phase as a result of thermo-chemi¬ 
cal reactions taking place in the chemically reactive bed 
of downdraft (biomass) gasifier. The local mass balance 
for each control volume in the fuel bed can be written as 

(Am,), = (m s ),_ x - (; m s ), = (m g ), - (m g ),_ x (3) 

where solid phase (s) comprises of biomass, char and ash, 
while gaseous phase ( g ) comprises of water vapour, vola¬ 
tiles, air and product gas. The solid conversion Am, in dry¬ 
ing and preheating zone is due to evaporation, thus, 
fraction of dry biomass, 7 db , and moisture content, F w , 
can be represented as 


Tdb + T w = 1 (4) 

In the pyrolysis zone, the mass fractions Y of volatile, char 
and ash can be related as 

*Z + 7 char + YZ = 1 (5) 

here, Y vyT denotes the mass fractions in dry wood and its 
subscripts viz., vol, char and ash denote the volatile, char 
and ash respectively. The solid conversion during the pyro¬ 
lysis process is obtained from the literature (Ragland and 
Aerts, 1991). During oxidation and reduction however 
the total char is consumed to form gases. Thus, 

yoxd , yered _ i (fit 

1 char T 1 char — 1 W 

here, and T^tr are the fractions of char consumed in 
the oxidation and reduction zones respectively. 

3.2.2. Momentum balance 

While the Darcy’s law for flow through different types of 
porous media is valid for Re 0p ^ 1, when Re Dp > 1, the 
inertial effects needs to be accounted for in the momentum 
balance. Thus, the correlations of Ergun (1952) and Mac¬ 
donald (Kaviany, 1995) as given in Table 1 are used in 
the present work for comparative analysis of the pressure 
drop through granular bed of gasifier. 

In Table 1, A' perm is the permeability of the medium, / is 
the thickness of bed, v 0 (= m/(Ap)) is the superficial veloc¬ 
ity and /h'| )p (= pv 0 d p / //) is the Reynolds number for par¬ 
ticle. The general equation to predict the pressure drop 
through the 7th CV in gasifier bed in terms of average fluid 


Table 1 

Correlations for prediction of pressure drop in porous bed 

S. no. Correlations Reference 


I 


II 


A P 

~—r= 


(1 - Sbed) 


+ hp4~. 3 - 

^bed 


^bed) 

d p 


A P 

~T 


= «MP!’o 


(1 — Sbed) 
Re D P d blsi dp 


+ PmP v I- 


Ergun (1952) [a = 150; /} = 1.75] 


Macdonald (Kaviany, 1995) [a = 180; /}= 1.8-4.0] 
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flow rate (w gm ), particle size and fluid density may be 
expressed as 


A Pi = 


( «(1 - £bed fli{Ti)li 



/77gm 


+ 


^bed )^7 



ra 


2 

gm 


(V) 


flow, Blasius’ relation gives, / fd = 0.316 /Re 015 for 
Re < 10 s (Fox and McDonald, 1995). The developing 
length for laminar flow is given by L Dev = 0.06 Re Z) Tuy , 
while for turbulent flow Z Dev = 4AReD l T ^ y . Average pres¬ 
sure drop parameter from entrance to asymptotic region 
of fully developed flow denoted by A" nev is given by the fol¬ 
lowing relation obtained by curve fitting to the data of 
Schmidt and Zeldin (White, 1991) for laminar pipe flow 


here, a and [’> are the two constants the values of which are 
given in Table 1. Macdonald fixed /? at 1.8 for smooth sur¬ 
face while 4 for rough surface. '/' s is the particles sphericity, 
A m is the average cross-sectional area of gasifier tube in 
each CV. In Eq. (7), the Mehta and Flawely (Churchill, 
1988) proposed the correction factor £ to compensate the 
drag of the confining wall by introducing surface area of 
the wall. 

f=[l + (2/3)(rf p /A)/(l-ebed)] (7a) 


K Dev = exp ( 0.3 — 2.9 x 10 j 

( ^Dev 


^Dev 


- 1.43 x 10~ 


\ReD 


Tuy. 


R.CD'Yuy 

2 \ 


(9) 


The Darcy friction factor for the developing region is then 
modified as 


/Tuy — /fd + 


K 


Dev 


i-Tuy/Dxuy 


( 10 ) 


Around the periphery of oxidation zone of the IISc-NET- 
PRO gasifier, the three tuyers are radially distributed at an 
angle of 120° to supply the air for gasification as shown in 
Fig. 3. These tuyers are made of straight pipe. The pressure 
drop across a tuyer is modeled as the developing regime for 
nearly the range of under consideration as given by 


AP-Tuy — 


A/uy +/ti lyE-ruy/Thuy 
v 2/r l]r /i Tuy y 


‘Tuy 


( 8 ) 


here, K Tuy is the pressure loss parameter to represent the 
combined effect of entrance and exit for the tuyer. 

In laminar fully developed flow (Re <2300), Darcy- 
Weisbach friction-factor, / fd is 64 /Re while for turbulent 



Fig. 3. Top view of the gasifier. 


3.2.3. Energy balance 

The axial and radial heat transport in gasifier bed 
accounted for advective, conductive and radiative heat 
fluxes at the boundaries of each control volume in the por¬ 
ous bed. The axial heat transfer accounted for thermal 
interactions between the neighboring CVs, while the radial 
heat loss is describing the heat loss to the surroundings. 
The advection is due to the flow of gas as well as bio¬ 
mass/char down the bed of the gasifier as depicted in 
Fig. 4. Thus, for each CV of each zone, the energy equation 
is written as 


[m s ll s + /77g/7 g ] / _ | + Q di ff 7 _i + Q gen j QlossJ 

= [m s h s + m g h g \j + Q dm (11) 


here h = f CpdT, is the sensible enthalpy change with re¬ 
spect to surroundings, it comprises of solid phase (biomass, 
char and ash) and gaseous phase (water vapour, volatiles, 
air and combustion products). 

Heat loss to the surroundings, <2i oss , is modeled in terms 
of equivalent thermal resistance, R s , as 


T, - T A _ Tt - T A 

Rs.t R ?(ins,7) H - -^/(air,/) 


( 12 ) 


Here R r{ ins ,z) and Rkaw.d are the resistance due to insulation 
at wall and outer surface of 7th CV in the bed. Qnsj-\ and 
Qdiff.i denote the heat blast in and out, and are modeled as 


Insulated reactor wall 


. 

dLj.i 

, 


CVtl (m, h s )C 

1 Rb,.i 

w 

§> ( m s hgh-i 

n 1 

i 

/ 

dU 


CVl Qgen.I T[( 

| Rbi 

\ Y ’ R., 

<)- 

0 | 

i 

cnv® Ta 

L> Qloss, 1 

dLi+i 


CV I+1 (m]h s ), i 

9 Qdifi (m g h g )j 

i+i 


Fig. 4. Energy balance for the control volume. 
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6diff,/-l — ( KsAVT)^ 


QdiffJ — 


T, - Tj +1 

Rb,I 


T ,-1 - T, 


(13) 


/?b,r and i?b,/-i are the thermal resistances in axial direction 
due to thermal interaction of CV/ with upper CV/_i and 
lower CV/+i. denotes the heat generation and depends 
on the process in a control volume: drying results in absorp¬ 
tion of latent heat of vaporization, pyrolysis results in heat 
of pyrolysis, oxidation or reduction processes results the 
heat generation or absorption due to gasification process. 


3.2.4. Effective thermal conductivity 

The sub-model of effective thermal conductivity, sche¬ 
matically represented in Fig. 5 has been described here 
for clarity. 

The contribution of conduction transport in the effective 
thermal conductivity follows the model of Sharma (2005), 
while contribution of radiation transport k rild has been 
adapted from the model of Sharma (2006). Here a cylindri¬ 
cal unit cell consists of two half spherical particles contact¬ 
ing each other with 6 ct as the half angle of contact. Heat 
transfer is assumed to be in the vertical direction parallel 
to the axis connecting the two spheres. Two conductive 
paths are considered between the two - the central region 
of the width of the contact surface, and the surrounding 
region where the heat transfer passage is from solid particle 
to the fluid in the void and then again through the solid 
particle, while the third path describes the radiation trans¬ 
port between the bed particles. For the model formulation, 
the particles are assumed to be spherical with equal diam¬ 
eter. The overall temperature drop, AT, in the cylindrical 
unit cell given as 

AT= AT s + AT g + AT' s (14) 

The heat conduction through solid-fluid-solid region in 
the unit cell can be represented by 


here, k effsg represents the contribution of conduction in so¬ 
lid-fluid-solid region to the effective thermal conductivity; 
A c i is the cross-sectional area of the unit cell and l c \ is its 
length, with A d = nflp/4; l d = d p and r = (<7 p /2)Sin0. Inte¬ 
grating Eq. (15) for limit, 9 = 9 ct to 0 max = ti/ 2 and assum¬ 
ing d v is constant within a unit cell A' cll ; S g is given by 

* eff ' ,g_ * g (hiC 2 + C 1 )-Jfc s C 1 (16) 


where C\ = Cos 0 ct ; C 2 = 1 — Cos 0 cl : while 9 ct is based on 
particle contact following Enoeda (Slavin et al., 2002) 
who took the contact radii to be typically 0.7% of particle 
diameter. Similarly, determining the conduction through 
the contact region, k eBcon is given by 


k 


eff,con — 



(17) 


According to Sharma (2006), the contribution of radiation 
in the cylindrical unit cell can be modeled assuming the 
wavelength-independent properties and for AT< 200 K, 
Tein and Drolen (1978) described the radiant conductivity 
assuming one-dinrensional, plane geometry with emitting 
particles in steady state, a deep bed with zero reflectivity 
at wall. Thus 

16 err 3 

K ff, rad = 7 = 4 crQdpT* (18) 


where Q is the shape factor and can be calculated in terms 
of extinction coefficient, ( a ex ). Kamiuto and Yee (2005) 
proposed the scaled extinction coefficient, ( 077 ) for large 
opaque particles as 



7=1 + 1-5(1 - Sbed) - 


1.5(27-l)(l- £bed ) 

dp 


(19) 


0.75(1 - £ b e d ) 2 for £ bed < 0.921 


( 20 ) 


j' 1 2nrdr % 2nrdr 

Ts Je a {dp/2) cos 0 g Ts Jg a d p (l -cos 9) 


^eff.sg^cl 


A d AT 


l, 


cl 


(15) 


For optically thick packed beds the internal radiant con¬ 
ductivity can be treated as a diffusive process. Thus, the 
effective thermal conductivity, k eB , can be written as 


keff — k, 


eff,sg ‘ 


( 21 ) 



Conduction 

(solids) 

Contactive 

Conduction 

(fluid) 

Radiation 

(voids) 


Fig. 5. Model of two contacting solid particles in the unit cell. 
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3.2.5. Particle size variation 

Here the initial diameter of biomass ( d h ) is an input 
parameter. In drying and preheating zones, there is no 
change in the diameter. However in pyrolysis, oxidation 
and reduction zones, feedstock undergoes chemical reac¬ 
tions leading to change in the particle size. As wood con¬ 
verts into char due to pyrolysis, its size reduces and its 
diameter at exit of pyrolysis zone can be given by 

dZr = MYZrPjPc har ) V3 ( 22 ) 

In oxidation zone, char is consumed to form gases. Particle 
size of char at the exit of this zone is given by 

<£ = «(! -Car)]' 73 (23) 

where p h and p tbar are the densities of the dry biomass par¬ 
ticle and the char particle respectively. In the reduction 
zone, all char is consumed and only ash particles (1 mm 
diameter size) leave the reduction zone. 

3.2.6. Bed porosity variation 

Bulk porosity of bed, e be d is determined using the corre¬ 
lation suggested by Chen and Gunkel (1987) 

«W = 0.5-0.2^1-0 (24a) 

Here d p is the particle diameter in the current control 
volume. 

Winterberg and Tsotsas, (2000) suggested the correla¬ 
tion for variation in radial porosity for spherical particles 
bed as 

£bed(r) = £ b ed(l + Cj.e Ct ^ (24b) 

where C 3 = 0.65/£ bed — 1, C 4 = 6.0 

The averaged bed porosity £ bed can be obtained by 
integration 

2 f R 

^bed — ~2 I ^bed (r)rdr (24c) 

B Jo 

4. Solution procedure 

The producer gas flow rate, air-fuel ratio and heat gen¬ 
eration in oxidation zone are required as input to the 
model. Then biomass consumption rate, water vapour 
and air inflow from open top as well as through the tuyers 
can be related each other using the energy balance Eq. (11) 
and momentum balance Eq. (7) for pressure drop in bed 
and Eq. ( 8 ) for pressure drop through the tuyers are cou¬ 
pled with global mass balance Eq. (2) and local mass bal¬ 
ance Eq. (3). Thus the sum of pressure drop across 
preheating, drying and pyrolysis zones can be related to 
pressure drop across the tuyers as 

APpreheat T AP dr y -T APpy ro — APj U y (25) 

here subscripts preheat, dry and pyro denote the preheat¬ 
ing, drying and pyrolysis zone respectively. The set of these 


equations for each 13 control volumes in conjunction with 
Eq. (25) is solved using Gauss-Siedel iteration of conver¬ 
gence to predict bed temperature profile and pressure drop 
across the gasifier bed. The total pressure drop across the 
gasifier thus is obtained as 

A P g = AP preheat + AP dry + AP p yr o + AP oxd + AP red (26) 

Here subscripts oxd and red denote the oxidation and 
reduction zone respectively. 

5. Property data 

The viscosity of gaseous phase is obtained from the cor¬ 
relation of Hagge and Bryden (2002) as given by 

p(T) = AMI x 10- 7 P 064487 (kgnr 1 s~‘) (27) 

Ragland and Aerts (1991) reported the thermal conductiv¬ 
ity (W/m K) of char and biomass as given by 

£char = 0.675 , cha r - 0.071 

k b = S b (0.1941 + 0.4064T w ) + 0.1864 + 0.002(P - T A ) 

(28a) 

k g st k dh = 3.0 x 10~ 4 r°' 79 (28b) 

here, S is the specific gravity (kJ/kg K) and its subscripts 
denote the biomass and char. Their values are obtained 
from literature (Ragland and Aerts, 1991). The data from 
Keenan et al. (Borman and Ragland, 1998) has been used 
to curve fit for thermal conductivity and specific heat of 
air and same curve fits are considered for gaseous phase, 
the relation is valid for 300 < T < 1500 K. 

The specific heat for working fluid (air) is obtained by 
curve fit to data of Keenan (Borman and Ragland, 1998) 

Cp hd = 0.1031 +0.003867P 

C Ph = [ 1 C Pbd + 4.19T W ]/(1 + Y v ) + 0.02355 T 

- 1.327 w -6.27 w (29a) 

Cp ch . dr = 1.39 + 0.00036P 

Cp s st Q? air = 0.931 + 0.0002P - 1.01 x 10~ 8 P 2 (29b) 

However, the Cp g and k g values of product gas can be ob¬ 
tained from the gas composition. 

Since the gaseous oxidation is very fast as compared to 
char oxidation (Ragland and Aerts, 1991) thus, the con¬ 
sumption of char in oxidation zone can be neglected in 
comparison with the char consumed in reduction zone 

oxd ^ yred \ 
char 1 char/* 

6. Simulation results and discussion 

The above model predicts the temperature profile, pres¬ 
sure drop through bed, airflow rates and apportionment 
of air inflow from the open top and through the tuyers for 
given mass flow rate of producer gas, air-fuel ratio, feed¬ 
stock size and heat generation or absorption in each zone. 
Jayah et al. (2003) reported experimentally that the air-fuel 
ratio for gasification process does not vary much (1.9-2.3 
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for wide range of gas flow rates). Thus, in the present work, 
the air-fuel ratio is maintained at 2.3. A realistic model 
requires the complete kinetics and energetics of the pyroly¬ 
sis, oxidation and reduction zones to be solved. The scope 
of this paper does not include these models. In order to 
account for actual hot flow conditions; the heat of pyrolysis, 
heat generation in oxidation zone and heat absorbed in 
reduction zone of the downdraft biomass gasifier are sup¬ 
plied as input parameters to the fluid flow and heat transfer 
model. For the present work, however, the overall pyrolysis 
is considered to be thermally neutral (Ragland and Aerts, 
1991) and the ratio of Q oxd and Q red is fixed at 2.8 (Sheng, 
1989). Using these approximations, the only parameter Q oxd 
would be needed as input information to describe the overall 
chemistry of chemical reactions in the gasifier. 

6.1. Validation 


10 -i 


TO 

o 


9 - 
8 - 
7 - 
6 - 
5 - 
4 - 
3 - 
2 - 
1 - 
0 -■ 


3 


A Ergun Equation (dp=40mm) 
o Ergun Equation (dp=38mm) 
• Experiments 


8 


o 

ft 


5 7 9 

Gas flow rate (g/s) 


8 

A 


11 


Fig. 6. Comparing the predictions of gasifier pressure drop with 
measurements for various gas flow rate, T w = 0.11, cf> = 0.4, Ergun 
correlation. 


It is quite difficult to quantify the parameter Q oxd from 
the traditional experiments on the gasifier. Thus, in order 
to validate the above fluid flow and heat transport model, 
the temperature profile may be supplied as input informa¬ 
tion to compare the predictions for the pressure drop 
through the gasifier bed against experimental data. The 
experiments are performed to obtain the temperature pro¬ 
files and gasifier pressure drops through the gasifier bed at 
different gas flow rate (described in Section 2: System 
description and experiments). These temperature profiles 
are converted into zonal temperature specifications for var¬ 
ious flow rates as listed in Table 2. The height of preheat¬ 
ing, drying, pyrolysis, oxidation and reduction zone are 
taken as 100, 60, 7, 17 and 17 cm respectively. 

With these temperatures specifications of each zone at 
the corresponding gas flow rates (Table 2), zone heights, 
with varying particle size distribution (Eqs. (22) and 
(23)); the simulations are performed for typical particle size 
of 38 ± 2 mm. The model predictions for pressure drop 
through bed in hot flow conditions are compared with 
experimental data as shown in the Fig. 6. The good agree¬ 
ment between predictions for pressure drop and experi¬ 
ments for actual operating condition is obtained. As 
expected, the pressure drop is found to be increased with 
increasing gas flow rates. 

6.2. Parametric analysis 

After validation of the model, the parametric analysis 
has been carried out to investigate the effect of parameters 

Table 2 


Zonal temperature specifications obtained from experiments 


»'pg (g/s) 

•^preheat 

-^drying 

T 

1 pyr 

T>xd 

Aed 

10.6 

301 

303 

880 

1270 

945 

8.5 

301 

303 

870 

1207 

940 

6 

301 

303 

860 

1204 

930 

4.7 

301 

304 

840 

1150 

920 


on the temperature profile and pressure drop through the 
gasifier bed. The baseline values for heat generation in oxi¬ 
dation zone, Q oxd , gas flow rate at exit, m pg , and air-fuel 
ratio are fixed at ll.OkW, 8 g/s, 2.3, respectively. With 
these baseline values and for spatially varying particle size 
distribution for d b value of 3 cm, the simulations are per¬ 
formed to investigate the effects of parameters - effective 
thermal conductivity, gas flow rate, heat released in oxida¬ 
tion zone, moisture content, particle size, air inlet temper¬ 
ature and thickness of insulation at reactor wall on the 
temperature profile and pressure drop through the gasifier 
bed. 

6.2.1. Effective thermal conductivity 

The predictions of temperature profile using the effective 
thermal conductivity model described in the present work 
have been compared with Petersen et al. (1997) and Slavin 
et al. (2002) as shown in Fig. 7, while Fig. 8 shows the pre¬ 
dictions for gasifier pressure drop through gasifier bed. The 
graphs reveal that the case with higher temperature profile 
in the gasifier bed predicts slightly higher value of pressure 
drop. At higher temperature the density of gases is lower, 



Fig. 7. Effect of k e s on bed temperature profile, m pg = 8 g/s, 
2oxd = 11 kW, 4 = 3 cm. 
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keff (W/m-K) 

Fig. 8. Effect of k e n on the pressure drop, ra pg = 8 g/s, goxd = 11 kW, 
d b = 3 cm, Ergun correlation. 

resulting in slightly higher velocities for given gas flow rate. 
Increased velocity and viscosity due to rise in temperatures, 
results in relatively higher pressure drop. If the temperature 
is doubled in a given CV, the fluid density is halved result¬ 
ing in doubling of the fluid velocities, for a given mass flow 
rate of gas. Since there is a quadratic dependence of pres¬ 
sure drop on velocity (Table 1), the corresponding increase 
in pressure drop in a given control volume will be nearly 
fourfold. Similarly the dependence of gas viscosity follows 
power law with temperature as given by the Eq. (27). 

6.2.2. Gas flow rate 

The effect of variation in gas flow rate from 7 to 10 g/s 
has been studied on the bed temperature profile and pres¬ 
sure drop through the bed for given heat generation in oxi¬ 
dation is shown in Figs. 9 and 10. Figs. 9 shows that the 
temperature of pyrolysis, oxidation and reduction zones 
decreases with increasing the gas flow rate due to greater 
advection of heat away from these zones. As expected, 
the pressure drop will be higher at high flow rate, since it 
offers higher resistance to flow due to increased velocity 
through the pores. Predictions for pressure drop also show 
that Macdonald equation for porous bed gives slightly 



Fig. 9. Effect of gas flow rate at gasifier exit, m m , Q ox & = 11 kW, 
d b = 3 cm. 



Gas flow rate (g/s) 

Fig. 10. Effect of gas flow rate at gasifier exit, m pg , Q OX d=llkW, 
d b = 3 cm. 

higher values of pressure drop as compared to Ergun 
equation. 

6.2.3. Heat generation in oxidation zone 

The effect of variation in heat generation in oxidation 
zone from 10 to 12 kW has been plotted against the tem¬ 
perature and pressure drop scale as shown in Figs. 11 
and 12. Fig. 11 shows that the temperature of the bed 
increases significantly with increase in heat generation in 
oxidation zone. An increase of 54 °C in oxidation zone 
temperature has been observed for increasing the heat gen¬ 
eration in oxidation zone from 10 to 11 kW at typical gas 
flow rate of 8 g/s. As expected, increasing the heat genera¬ 
tion in oxidation zone results higher temperature of the 
neighboring zones (downstream reduction zone and/or 
upstream pyrolysis zone), while the temperatures in the 
drying and preheating zones do not show any significant 
effect. Fig. 12, shows that the pressure drop across the gas¬ 
ifier bed increases marginally with increase in heat genera¬ 
tion in oxidation zone. As the temperatures in oxidation 
and reduction zones improve significantly with heat gener¬ 
ation in oxidation zone, leading to increase in pressure 
drop through the gasifier bed. 



d b = 3 cm. 
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Fig. 12. Effect of heat generation in oxidation m pg = 8 g/s, d^ = 3 cm. 



- - Ergun Equation 

— Macdonald equation 


1 4- T - T -i— 

0 0.05 0.1 0.15 


Moisture content in biomass (drybasis) 


Fig. 14. Effect of moisture content in biomass, m pg = 8 g/s, Q oxd = 11 kW, 
db = 3 cm. 


6.2.4. Moisture content 

The moisture content in biomass feedstock is a critical 
operating parameter and the success of gasification opera¬ 
tion depends on it. The moisture content of sun dried wood 
used as biomass feedstock varies up to 20% dry basis (Rag¬ 
land and Aerts, 1991). The effect of variation of moisture 
content in dry biomass on the bed temperature and pres¬ 
sure drop is shown in Figs. 13 and 14. As moisture content 
in dry biomass increases, the temperature near the drying 
zone decreases (5). It could be expected due to the evapo¬ 
ration and superheating of the water vapour, while the tem¬ 
peratures of oxidation and reduction zone do not show any 
significant effect with variation in moisture content in feed¬ 
stock. Fig. 14 does not show any significant effect of mois¬ 
ture content on the pressure drop, since the variation of 
temperature in oxidation and reduction zone is insignifi¬ 
cant (5). 

6.2.5. Particle size 

The variation of initial feedstock size from 2 to 4 cm has 
been investigated on bed temperature profile and pressure 
drop through gasifier bed for spatially varying particle size 
distribution as shown in Figs. 15 and 16. Fig. 15 shows that 
the temperatures in oxidation and reduction zone improve 



Fig. 13. Effect of moisture content in biomass m pg = 8 g/s, g oxl j = 11 kW, 
d h = 3 cm. 
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Fig. 15. Effect of feedstock diameter, d b , on m pg = 8 g/s, Q oxd = 11 kW. 



Particle diameter (cm) 


Fig. 16. Effect of feedstock diameter, d b , on m pg = 8 g/s, 2oxd = 11 kW. 

significantly with reduction in feedstock size. The pressure 
drop in the bed decreases considerably with increase in par¬ 
ticle size (16). The effect of variation in particle size on tem¬ 
perature is understandable, since the contribution of 
radiation in the effective thermal conductivity of the bed 
increases with increase in particle size (Eq. (18)), thus, the 
heat blast from the oxidation zone to the neighboring zones 
(i.e. pyrolysis and reduction zone) increases as particle size 
increases since the heat generation in oxidation zone is 
independent of particle size, therefore, the temperature of 
oxidation zone start decreasing with increase in particle 
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size. The pressure drop in the gasifier bed also decreases, 
since the pore size increases with particle size and thus 
offers less resistance to flow with increase in particle size 
of feedstock. 

6.2.6. Input air temperature 

Jayah et al. (2003) demonstrated that the hot air 
improves the gasifier performance, as it affects the air-fuel 
ratio for the gasification operation. The effect of preheating 
the input air from 300 to 400 K has been investigated here 
as shown in Figs. 17 and 18. Fig. 17 shows that the temper¬ 
ature profile improves considerably with rise in tempera¬ 
ture of input air. The maximum temperature in the bed 
improves from 1133 to 1211 K for typical increase in air 
input temperature from 300 K to 400 K. Figs. 18 shows 
that pressure drop across the gasifier increases marginally 
with increase in input air temperature. Although the pres¬ 
ent gasifier design has been installed an annular jacket zone 
to recover the heat from the producer gas leaving the gas¬ 
ifier, yet the gas temperature even then is considerably 
high. The graph reveals that the length of heat recovery 
zone/annular jacket zone may be increased up to the pyro¬ 
lysis zone in order to enhance the preheating of input air. 



2oxd = 11 kw, rf b = 3 cm. 



Fig. 18. Effect of ambient temperature, T a, on the m pg = 8 g/s, 
2oxd = 11 kW, rf b = 3 cm. 


6.2. 7. Thickness of thermal insulation 

The lack of insulation thickness is the greatest contribu¬ 
tor to heat loss, which also determines the thermal effi¬ 
ciency of the reactor. The effect on temperature profile 
and pressure drop in the gasifier bed is plotted in Figs. 19 
and 20. As expected the temperature profile improves con¬ 
siderably with the initial increase in insulation thickness. 
An increase of about 47 °C in the temperature of oxidation 
and 73 °C in the reduction zone has been predicted for 
increasing the insulation of thickness at reactor wall from 
2 cm to 30 cm at the typical gas flow rate of 8 g/s. Fig. 20 
demonstrate the effect of insulation thickness on the pres¬ 
sure drop. The graph show the pressure drop increases 
sharply for initial increase in the insulation thickness there¬ 
after it become almost independent of insulation thickness. 
For the present design, the insulation thickness should not 
be less than 15 cm, otherwise it may affect the pressure 
drop and the reaction temperatures in the bed and hence 
the performance of the reactor. 



Fig. 19. Effect of insulation thickness on temperature profile 
2gen,4 = 11 kW, m pg = 8 g/s, 4 = 3 cm, Y„ = 0.11. 
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Fig. 20. Effect of insulation thickness on the pressure drop, 
2gen,4 = 11 kW, m pg = 8 g/s, F w = 0.11, 4 = 3 cm. 
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7. Conclusions 

Fluid flow and heat transport model for the reactive, 
porous bed of biomass gasifier has been presented for sim¬ 
ulating the apportionment of inflow air from open top and 
through tyurs, bed temperature profile, pressure drop 
through the gasifier bed. The conservation, fluid and heat 
transport model equations for biomass gasifier have been 
written for the six zones viz, preheating, annular jacket, 
drying, pyrolysis, oxidation and reduction zones; which 
comprises of total 13 CVs. Simulations are performed to 
predict the temperature profile and pressure drop for vary¬ 
ing various operating parameters, viz., gas flow rate, heat 
generation in oxidation zone, moisture content in feed¬ 
stock, particle size, and input air temperature. The simu¬ 
lated results show that the temperature profile is found to 
be very sensitive for gas flow rate, heat released in oxida¬ 
tion zone and input air temperature, while the temperature 
profile in the bed is found relatively insensitive for moisture 
content in feedstock. 

The results also show that the high gas flow and bed 
temperature; and the reduction in feedstock size is found 
to cause a marked increase in pressure drop through the 
gasifier, while the pressure drop is relatively insensitive to 
the moisture content. The simulation also revealed that 
the size of heat exchanger may be extended up to the pyro¬ 
lysis zone to enhance further preheating of the input air. 
The thickness of thermal insulation should not be less than 
15 cm at reactor wall. The simulation shows that 40-50% 
of the total air induced is inducted through the three tuyers 
and the rest comes in through the open top. This ratio does 
not vary much with variation in temperature, producer gas 
flow rate or particle size. 
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